The phase diagram of neutral quark matter: 
The effect of neutrino trapping 
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We study the effect of neutrino trapping on the phase diagram of dense, locally neutral three- flavor 
quark matter within the framework of a Nambu-Jona-Lasinio model. In the analysis, dynamically 
generated quark masses are taken into account self-consistently. The phase diagrams in the plane 
of temperature and quark chemical potential, as well as in the plane of temperature and lepton- 
number chemical potential are presented. We show that neutrino trapping favors two-flavor color 
superconductivity and disfavors the color-flavor-locked phase at intermediate densities of matter. 
At the same time, the location of the critical line separating the two-flavor color-superconducting 
phase and the normal phase of quark matter is little affected by the presence of neutrinos. The 
implications of these results for the evolution of protoneutron stars are briefly discussed. 
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I. INTRODUCTION 

Recently, interest in the structure of the phase diagram 
of neutral dense quark matter has increased considerably 
[H IE IE EI IE El ■ This was driven in part by several key 
observations regarding the influence of charge neutrality 
and equilibrium on the competition of various phases 
in quark matter 0, IE IE EE El El 

In general, dense quark matter is expected to be color- 
superconducting. (For reviews on color superconductiv- 
ity, see Ref. [l3|.) At extremely large densities, namely 
when the quark chemical potential /i is much larger 
than the constituent, medium-modified quark masses, the 
ground state of matter is give n by the so-called color- 
flavor-locked (CFL) phase fl4l ( for studies of QCD at 
large densities, see also Ref. [l5|). At the highest den- 
sities existing in stars, however, the chemical potential 
is unlikely to be much larger than 500 MeV, while the 
constituent mass of the strange quark is not smaller than 
the current mass, which is about 100 MeV. 

In stellar matter, therefore, the heavier strange quark 
may not be able to participate in diquark Cooper pair- 
ing as easily as the light up and down quarks. Then, the 
pairing of light quarks can lead to the two-flavor color- 
superconducting (2SC) ground state 0,0. It should be 
pointed out, however, that the 2SC phase is subject to 
large penalties after imposing the charge neutrality and 
the /3-equilibrium conditions 0] . Indeed, when the frac- 
tion of strange quarks and leptons is small, the electric 
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neutrality requires roughly twice as many down quarks 
as up quarks. In this case, Cooper pairing of up quarks 
with down quarks of opposite momenta becomes rather 
difficult. Then, depending on the details of the interac- 
tion, the gapless CFL (gCFL) phase E3i the gapless 2SC 
(g2SC) phase [Tflj . or even the normal quark (NQ) mat- 
ter phase may be more favored. Also, there exist other 
reasonable possibilities [18l lltA |2JJ, |2lj, |22j, |23| which, in 
view of the known instabilities in the gapless phases [24| , 
are considered to be very promising. 

In Ref. , we obtained the phase diagram of neutral 
quark matter, described by the NJL-type model with the 
parameter set from Ref. [25j . In contrast to previous stud- 
ies in Refs. 0, 0, 0, the dynamically generated quark 
masses were treated self-consistently in Ref. [j| . (For ear- 
lier studies on color superconductivity, treating quark 
masses as dynamical quantities, see also Refs. [2a. I^.) 
In the present article we follow the same approach to 
study the effect of a nonzero neutrino (or, more precisely, 
lepton-number) chemical potential on the structure of the 
phase diagram. This is expected to have a potential rele- 
vance for the physics of protoneutron stars where neutri- 
nos are trapped during the first few seconds of the stellar 
evolution. 

The effect of neutrino trapping on color- 
superconducting quark matter has been previously 
discussed in Ref. It was found that a nonzero 
neutrino chemical potential favors the 2SC phase and 
disfavors the CFL phase. This is not unexpected be- 
cause the neutrino chemical potential is related to the 
conserved lepton number in the system and therefore it 
also favors the presence of (negatively) charged leptons. 
This helps 2SC-type pairing because electrical neutrality 
in quark matter can be achieved without inducing a 
very large mismatch between the Fermi surfaces of up 
and down quarks. The CFL phase, on the other hand, is 
electrically and color neutral in the absence of charged 
leptons when T = [28| . A nonzero neutrino chemical 
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potential can only spoil CFL-type pairing. 

In the present paper we extend the analysis of Ref. Q 
by performing a more systematic survey of the phase di- 
agram in the space of temperature, quark and lepton- 
number chemical potentials. This also includes the pos- 
sibility of gapless phases, which have not been taken 
into account in Ref. @- While such phases are gener- 
ally unstable at zero temperature l24l. this is not always 
the case at nonzero temperature |29j. Keeping this in 
mind, we shall merely localize the "problematic" regions 
in the phase diagram, where unconventional pairing is 
unavoidable. We shall refrain, however, from speculating 
on various pos sibilities for the true ground state (see, e.g., 
since these are still under debate. 

In the application to protoneutron stars, it is of inter- 
est to cover a range of parameters that could provide a 
total lepton fraction in quark matter of up to about 0.4. 
This is the value of the lepton-to-baryon charge ratio in 
iron cores of progenitor stars. Because of the conserva- 
tion of both lepton and baryon charges, this value is also 
close to the lepton fraction in protoneutron stars at early 
times, when the leptons did not have a chance to dif- 
fuse through dense matter and escape from the star [23| • 
The densities in the central regions of protoneutron stars 
can reach values of up to about 10 times lar ger t han the 
nuclear saturation density po « 0.15 fm -3 |30| . In the 
model studied in the present paper, almost the whole 
range of possibilities shall be covered by restricting the 
quark chemical potential to fi < 500 MeV and the neu- 
trino chemical potential to fj, Ve < 400 MeV. 

This paper is organized as follows. In SeclHl we briefly 
introduce the model. The general effects of neutrino trap- 
ping on two- and three-flavor quark matter are discussed 
in Sec. IIIII The results for the phase diagram of quark 
matter with neutrino trapping are presented in Sec. IIVI 
In Sec. [Vj we summarize our results and discuss their 
implications for the evolution of protoneutron stars. 



grangian density reads 

8 

C = ip m)ip + G s ^2 [(''Ma'i/') 2 + (^75A a -0) 2 

a=0 

+ g d J2 [r^^^aU^] [(fc);*T5£ m £r SC C] 

- K jdet [0 (1 + 75 ) V] + det [4> (1 - 75 ) ^] J , (1) 

where the quark spinor field ip^ carries color (a = r,g,b) 
and flavor (a — u,<i, s) indices, rh = diagj(m„, m^, m s ) 
is the matrix of quark current masses; \ a (a = 1, . . . , 8) 
are the Gell-Mann matrices in flavor space, and Ao = 
^/2/3 1/. The charge conjugate spinor is defined as ipe — 
Cip T and -0c = ip T C, where -0 = ip^j is the Dirac con- 
jugate spinor and C = ij 2 j° is the charge conjugation 
matrix. 

In the model at hand, there are six mutually commut- 
ing conserved charge densities. They split naturally into 
the following three classes. 

(i) Electric charge. This is related to the U(l) symme- 
try of electromagnetism. The corresponding charge 
density is given by 

™Q = (0 f W) ~n e - , (2) 

where Q — diagy (|, — i, — -i) is the electric charge 
matrix of the quarks, and n e and n M denote the 
number densities of electrons and muons, respec- 
tively. 

(ii) Two lepton charges. As long as the neutrinos are 
trapped and their oscillations are neglected, the 
lepton family numbers are conserved. The corre- 
sponding densities read 

«L C = n e + n Ve , n Lft = n M + n VfJi . (3) 



II. MODEL 

We consider a system of up, down, and strange quarks, 
in weak equilibrium with charged leptons and the corre- 
sponding neutrinos. We assume that the lepton sector of 
the model is given by an ideal gas of massive electrons 
(m e ~ 0.511 MeV) and muons (m^ « 105.66 MeV), as 
well as massless electron and muon neutrinos. We do not 
take into account the r lepton, which is too heavy to play 
any role in dense matter. Also, we neglect the possibil- 
ity of neutrino mixing, and therefore we do not take into 
account the r neutrino either. 

In the quark sector, we use the same three-flavor 
Nambu-Jona-Lasinio model as in Ref. la. The La- 



(iii) Baryon number and two color charges. The SU (3)- 
color symmetry and U(l) baryon number symme- 
try imply the conservation of three independent 
charge densities in the quark sector, 

n=(il>H), n 3 = (^T 3 ij), n 8 = (^T 8 ^), (4) 

where T 3 = diag c (|, -\, 0) and V3r 8 = 
diag c (i, i, — 1) are the matrices associated with 
the two mutually commuting color charges of the 
SU (3) c gauge group. Note that n is the quark num- 
ber density that is related to the baryon number 
density as follows: ub = \n. An alternative choice 
of the three conserved charges is given by the num- 
ber densities of red, green, and blue quarks, i.e., 
n r = nB + n^ + ns/Vi, n g — ns — n 3 + n 8 /V3, and 
nb = n B - 2n 8 /\/3. 
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The six conserved charge densities defined above are re- 
lated to the six chemical potentials of the model. These 
are the quark chemical potential /i = /i_b/3, the two color 
chemical potentials 113 and /is, the electric charge chemi- 
cal potential /xq, and the two lepton-number chemical 
potentials and fJ,L ■ 

In chemical equilibrium, the chemical potentials of all 
individual quark and lepton species can be expressed in 
terms of these six chemical potentials according to their 
content of conserved charges. For the quarks, which carry 
quark number, color, and electric charge, this is related 
to the diagonal elements of the following matrix: 

1$ = + MQ Q^) <U + [M3 (T 3 ) ab + Ms (T 8 )J S?f> , 

(5) 

i.e., the chemical potential of a quark with flavor a and 
color a is given by ^" 
hand, carry lepton number only. Then 



The neutrinos, on the other 



(G) 



Finally, electrons and muons carry both, lepton number 
and electric charge, hence 



/V = Mi,, - 



(7) 



As in Ref. [£j , the quark part of our model is treated in 
the mean-field (Hartree) approximation, allowing for the 
presence of both quark-antiquark condensates and scalar 
diquark condensates. In the model, the constituent quark 
masses are defined by the following relation: 



4G s <r a + 2Ka /3 a 1 



(8) 



where a a = (ipa^a) ( n0 sum over a) are the quark- 
antiquark condensates, and the set of indices (a, (3, 7) 
is a permutation of (u, d, s). 

Diquark pairing in the model is described by three gap 
parameters 



2G D ($ c ) a a vy 5 e 



af3c 



(9) 



(no sum over c). By definition, the parameters Ai, A2, 
and A3 correspond to the down-strange, up-strange, and 
up-down diquark condensates, respectively. 

By following the same steps in the derivation as in 
Ref. [j| and including the ideal-gas contributions for the 
leptons, we arrive at the following expression for the pres- 
sure: 



1 

2^ 



E 



dk k 2 



kl + 2Tln 1 



(l + e-¥) 



4Ka u a d a s - — £ | A c | 2 - 2G S £ a 2 a 



E 

l—e,fi 

1 

24tt 2 



T 
Z2 



dk k 2 In 1 + e 



-) 



2n 2 ^T 2 + —n 4 T 



(10) 



where ej are eighteen independent positive-energy eigen- 
values, see Ref. 5J for details. The pressure in Eq. l(TU|) 
has a physical meaning only at stationary points 
with respect to variations of the chiral and color- 
superconducting condensates. Thus, a a and A c must sat- 
isfy the following set of six gap equations: 



da a 



dp 
8A C 







(11a) 
(lib) 



To enforce the conditions of local charge neutrality in 
dense matter, one also requires that the following three 
equations are satisfied: 



»3 



n 8 = 



dp 
dp 
dp 







= 



= 



(12a) 
(12b) 
(12c) 



We should note that, in general, the restriction to the 
color charges 713 and ns is not sufficient, but one has to 
ensure that the densities n a = (ip^T a i/)) vanish for all 
generators T a of SU(3) C , i.e., a = 1, ... ,8 [31|. However, 
for all phases we are considering here, we have checked 
explicitly that the condition n a — is automatically sat- 
isfied for o ^ 3,8. Thus, we are left with Eqs. l(rl|l. By 
solving these, we determine the values of the three corre- 
sponding chemical potentials [Aq, ^3, and ^8 for a given 
set of the other chemical potentials, fj,, (J-l c , fJ-L , and for 
a given temperature T. In general, therefore, the phase 
diagram of dense quark matter with neutrino trapping 
should span a four-dimensional parameter space. 

Note that instead of using the chemical potentials, fi, 
/j,L b , and jj,Lu. as free parameters in the study of the phase 
diagram, one may also try to utilize the quark number 
density and the two lepton fractions, 



Yl. 
Yl u 



dp 
dpi 

3^ 
n 



where the two lepton densities are defined by 



dp 
dfiL e 



dp 



(13a) 
(13b) 
(13c) 

(14) 



In some cases, the choice of n, Yl b and Yl as free para- 
meters is indeed very useful. For instance, this is helpful 
in order to determine the initial state of matter inside the 
protoneutron star at very early times, when the lepton 
fractions are approximately the same as in the progenitor 
star (i.e., Yj, e « 0.4 and YV = 0). The problem is, how- 
ever, that such an approach becomes ambiguous in the 
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vicinity of first-order phase transitions, where the baryon 
number density as well as the lepton fractions are in gene- 
ral discontinuous. For this reason, it is more appropriate 
to study the phase structure of (dense) QCD at given 
fixed values of the chemical potentials fj,, Hl b i and //£ . 
Unlike densities, the chemical potentials change continu- 
ously when the system crosses a boundary of a first-order 
phase transition. (Here it should be noted that, because 
of the long-range Coulomb interaction enforcing the con- 
straints tiq = 0, n 3 = and n 8 = 0, the chemical poten- 
tials fiQ, /^3, and may change almost discontinuously 
at a boundary of a first-order phase transition.) 

III. SIMPLIFIED CONSIDERATIONS 

As mentioned in the Introduction, neutrino trapping 
favors the 2SC phase and strongly disfavors the CFL 
phase 0. This is a consequence of the modified (3- 
equilibrium condition in the system. In this section, we 
would like to emphasize that this is a model-independent 
effect. In order to understand the physics behind it, it is 
instructive to start our consideration from a very simple 
toy model. Later, many of its qualitative features will be 
also observed in our self-consistent numerical analysis of 
the NJL model. 

Let us first assume that strange quarks are very heavy 
and consider a gas of non-interacting massless up and 
down quarks in the normal phase at T — 0. As required 
by (3 equilibrium, electrons and electron neutrinos are 
also present in the system. (Note that in this section we 
neglect muons and muon neutrinos for simplicity.) 

In the absence of Cooper pairing, the densities of 
quarks and leptons are given by 

n u ,d = — o- , n e = r-5 , n Ve = . (15) 

Expressing the chemical potentials through /i, fj,Q and 
/iL e , and imposing electric charge neutrality, one arrives 
at the following relation: 

2(1 + \vf (1 - \yf (x - yf = , (16) 

where we have introduced the chemical potential ratios 
x = /UL e //i and y = /j,q/(i. The above cubic equation can 
be solved for y (electric chemical potential) at any given 
x (lepton- number chemical potential). The result can be 
used to calculate the ratio of quark chemical potentials, 

Md/ Mu = ( 3 - y)/( 3 + 2 v)- 

The ratio [id/fJ-u as a function of fi^^/fi is shown in 
Fig.ffl At vanishing one finds y w —0.219 and, thus, 
Hd/ A*« ~ 1.256 (note that this value is very close to 2 1 / 3 w 
1.260). This result corresponds to the following ratios of 
the number densities in the system: n u /n,d ~ 0.504 and 
n e /nd « 0.003, reflecting that the density of electrons is 
tiny and the charge of the up quarks has to be balanced 
by approximately twice as many down quarks. 



At = /i, on the other hand, the real solution to 
Eq. (|16fl is y = 0, i.e., the up and down Fermi momenta 
become equal. This can be seen most easily if one inverts 
the problem and solves Eq. I|16l) for x at given y. When 
y = one finds x = 1, meaning that fid = \i u and, in turn, 
suggesting that pairing between up and down quarks is 
unobstructed at /i£ e = fi. This is in contrast to the case 
of vanishing fiL e , when the two Fermi surfaces are split 
by about 25%, and pairing is difficult. 

It is appropriate to mention that many features of the 
above considerations would not change much even when 
Cooper pairing is taken into account. The reason is that 
the corresponding corrections to the quark densities are 
parametrically suppressed by a factor of order (A/ju) 2 . 

In order to estimate the magnitude of the effect in 
the case of dense matter in protoneutron stars, we in- 
dicate several typical values of the lepton fractions Yl c 
in Fig. ^ As mentioned earlier, Yl c is expected to be 
of order 0.4 right after the collapse of the iron core of 
a progenitor star. According to Fig. ^ this corresponds 
to fid/fJ-u ~ 1-1, i-C, while the splitting between the up 
and down Fermi surfaces does not disappear completely, 
it gets reduced considerably compared to its value in the 
absence of trapped neutrinos. This reduction substan- 
tially facilitates the cross-flavor pairing of up and down 
quarks. The effect is gradually washed out during about 
a dozen of seconds of the deleptonization period when 
the value of Yj, e decreases to zero. 

The toy model is easily modified to the opposite ex- 
treme of three massless quark flavors. Basically, this cor- 
responds to replacing Eq. (|16l) by 

2(1 + ^) 3 - 2(1 -iy) 3 - [x-yf = . (17) 

In the absence of neutrino trapping, x = 0, the only real 
solution to this equation is y = 0, indicating that the 
chemical potentials (which also coincide with the Fermi 
momenta) of up, down, and strange quarks are equal. 
This reflects the fact that the system with equal densi- 
ties of up, down, and strange quarks is neutral by itself, 
without electrons. With increasing x oc /il c , the solution 
requires a nonzero y oc /iq, suggesting that up-down and 
up-strange pairing becomes more difficult. To see this 
more clearly, we can go one step further in the analysis 
of the toy model. 

Let us assume that the quarks are paired in a regular, 
i.e., fully gapped, CFL phase at T = 0. Then, as shown 
in Ref. J28J, the quark part of the matter is automati- 
cally electrically neutral. Hence, if we want to keep the 
whole system electrically and color neutral, there must 
be no electrons. Obviously, this is easily realized with- 
out trapped neutrinos by setting /xq equal to zero. At 
non- vanishing fi^ c the situation is more complicated. The 
quark part is still neutral by itself and therefore no elec- 
trons are admitted. Hence, the electron chemical poten- 
tial n e = hl c — must vanish, and consequently /iq 
should be nonzero and equal to \ih B - It is natural to ask 
what should be the values of the color chemical potentials 
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/X3 and (is in the CFL phase when //L e 7^ 0. 

In order to analyze the stress on the CFL phase due to 
nonzero /i£ e , we follow the same approach as in Refs. 0, 
ITlj . In this analytical consideration, we also account for 
the effect of the strange quark mass simply by shifting 
the strange quark chemical potential by — Af 2 /(2/i). In 
our notation, pairing of CFL-type requires the following 
"common" values of the Fermi momenta of paired quarks: 



common 
PF,(ru,gd,bs) 

common 
l J F,(rd,gu) 

common 
PF,(rs,bu) 

common 
PF,(gs,bd) 



M - 

M ' 
M' 
/' 



Ml 
mq 

6 

Mq 
6 

MQ 
3 



Ms 

2^3' 
4 



M3 

4 



Ms 

Ms 
4^3 



(18a) 
(18b) 

M 2 

(18d) 



These are used to calculate the pressure in the toy model, 



.(toy) 



^ 3 3 
^2 X/ X/ 



a—1 a—1 
,2 A2 



Pf,c 



+3 



(Ma ~P)p 2 dp 
(ml c ~ Mq) 4 . Mi e 



12tt 2 



24tt 2 ' 



(19) 



By making use of this expression, one easily derives the 
neutrality conditions as in Eq. (|12|) . In order to solve 
them, it is useful to note that 



«3 



1 

7! 



(MQ ~ ML e 
3tt 2 



(20) 



Thus, it becomes obvious that charge neutrality requires 
MQ — Mi e - The other useful observation is that the ex- 
pression for n 3 is proportional to ^3 + /zq. So, it is van- 
ishing if (and only if) /i 3 = —fJ-Q, which means that 
/J3 = — ML e - Finally, one can check that the third neu- 
trality condition n s = requires 



Ms 



Air 



\/3 "v/3/i 



(21) 



The results for the charge chemical potentials [Iq, /i 3 , and 
fi$ imply the following magnitude of stress on pairing in 
the CFL phase: 



3^(gs,bd) 







Mr 




2 




M& 




Mr 




2 










2 



A/ 2 

= Ml 



(22a) 
(22b) 

(22c) 



Note that there is no mismatch between the values of 
the chemical potentials of the other three quarks, /i" = 
/^ = Mo=M~M s 7(6m)- 



From Eq. (|22|l we see that the largest mismatch occurs 
in the (rs,bu) pair (for positive M£ e )- The CFL phase 
can withstand the stress only if the value of Sfi^ rs j, u ^ is 
less than A 2 . A larger mismatch should drive a transition 
to a gapless phase exactly as in Refs. 0,0]. Thus, the 
critical value of the lepton-number chemical potential is 



(cr) 



Ml 

2fi ' 



(23) 



When Mi e > Ml^? the CFL phase turns into the gCFL' 
phase, which is a variant of the gCFL phase ^l| . By defi- 
nition, the gapless mode with a linear dispersion relation 
in the gCFL' phase is rs-bu instead of gs-bd as in the 
standard gCFL phase. (Let us remind that the mode aa- 
bfl is defined by its dispersion relation which interpolates 
between the dispersion relations of hole-type excitations 
of aa-quark at small momenta, k <C /1" , and particle-type 
excitations of 6/3-quark at large momenta, k 3> Mfc ■) 

In order to see what this means for the physics of pro- 
toneutron stars, we should again try to relate the value of 
/i£ e to the lepton fraction. As we have seen, there are no 
electrons in the (regular) CFL phase at T = 0. Therefore 
the entire lepton number is carried by neutrinos. For the 
baryon density we may neglect the pairing effects to first 
approximation and employ the ideal-gas relations. This 
yields 



Yr 



6 ^ fi 



(24) 



Inserting typical numbers, /1 = 500 MeV and < A « 
50-100 MeV, one finds Y Le < lO^-lO" 3 . Thus' there is 
practically no chance to find a sizable amount of leptons 
in the CFL phase. The constraint gets relaxed slightly 
at nonzero temperatures and/or in the gCFL phase, but 
the lepton fraction remains rather small even then (our 
numerical results indicate that, in general, Y^ e < 0.05 in 
the CFL phase). 



IV. RESULTS 

The simple toy-model considerations in the previous 
section give a qualitative understanding of the effect of 
neutrino trapping on the mismatch of the quark Fermi 
momenta and, thus, on the pairing properties of two- 
and three-flavor quark matter. Now, we turn to a more 
detailed numerical analysis of the phase diagram in the 
NJL model defined in Sec. [H] 

In the numerical calculations, we use the following set 
of model parameters [2f| : 



A = 602.3 MeV , 
GsA 2 = 1.835 , 
KA 5 = 12.36 , 

m s = 140.7 MeV , 
m ud = 5.5 MeV . 



(25a) 
(25b) 
(25c) 
(25d) 
(25e) 
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The parameters are chosen to reproduce several key ob- 
servables of vacuum QCD [25|. In this paper, our choice 
for the diquark coupling is Gd — jGs, corresponding to 
what is called "intermediate coupling" in Ref. 

In order to obtain the phase diagram, one has to de- 
termine the ground state of matter for each given set 
of the parameters. As discussed in Sec. [n] in the case 
of locally neutral matter with trapped neutrinos, there 
are four parameters that should be specified: the tempe- 
rature T, the quark chemical potential fj, as well as the 
two lepton family chemical potentials and . After 
these are fixed, the values of the pressure in all competing 
neutral phases of quark matter should be compared. This 
is determined by using the same algorithm as in Ref. p| . 
The complete set of equations (1111) and ljl2|l is solved for 
each of the eight phases allowed by symmetries. Then, 
the corresponding values of the pressure are determined 
from Eq. (|l(Jf) . The phase with the largest pressure is the 
ground state. 

In the present paper, we always assume that the muon 
lcpton-number chemical potential vanishes, i.e., /i£ fi = 0. 
This is expected to be a good approximation for mat- 
ter inside a protoneutron star. Our analysis can thus 
be interpreted as an extension of the T-[i phase dia- 
gram discussed in Ref. [f| into the direction. Con- 
sequently, the complete phase structure requires a three- 
dimensional presentation. 



A. Three-dimensional phase diagram 

The general features of the phase diagram in the three- 
dimensional space spanned by the quark chemical poten- 
tial /it, the lepton- number chemical potential /i£ c , and the 
temperature T are depicted in Fig. [21 Because of a rather 
complicated structure of the diagram, only the four main 
phases (xSB, NQ, 2SC and CFL) are shown explicitly. 
Although it is not labeled, a thin slice of a fifth phase, 
the uSC phase, squeezed in between the 2SC and CFL 
phases, can also be seen. While lacking detailed informa- 
tion, the phase diagram in Fig. gives a clear overall 
picture. Among other things, one sees, for example, that 
the CFL phase becomes strongly disfavored with increas- 
ing fj,L c and gets gradually replaced by the 2SC phase. 

In order to discuss the structure of the phase dia- 
gram in more detail we proceed by showing several two- 
dimensional slices. These are obtained by keeping one of 
the chemical potentials, /i or fixed and varying the 
other two parameters. 



B. T—/J, phase diagram 

We begin with presenting the phase diagrams at three 
fixed values of the lepton-number chemical potential, 
/ii e = (upper panel), /*L e = 200 MeV (middle panel) 
and /i Le = 400 MeV (lower panel), in Fig. The first 
of them has already been discussed in Ref. [5| . Here we 



show it again only for the purpose of comparison with 
the other two diagrams. The general effect of neutrino 
trapping can be understood by analyzing the similarities 
and differences between the three diagrams. Note that, 
in this paper, we use the same convention for line styles 
as in Ref. [f|: thick and thin solid lines denote first- and 
second-order phase transitions, respectively; dashed lines 
indicate the (dis-)appearance of gapless modes in diffe- 
rent phases. 

Here it is appropriate to note that, in the same model, 
a schematic version of the T-/i phase diagram at /i£ c = 
200 MeV was first presented in Ref. [g, see the right 
panel of Fig. 4 there. If one ignores the complications 
due to the presence of the uSC phase and various gapless 
phases, the results of Ref. || are in qualitative agreement 
with our findings here. 

In order to understand the basic characteristics of 
different phases in the phase diagrams in Fig. [3J we 
also present the results for the dynamical quark masses, 
the gap parameters, and the charge chemical potentials. 
These are plotted as functions of the quark chemical 
potential in Figs. 0] and [5J for two different values of 
the temperature in the case of /i£ e = 200 MeV and 
/ii e = 400 MeV, respectively. 

In each of the diagrams, there are roughly four dis- 
tinct regimes. At low temperature and low quark che- 
mical potential, there is a region in which the approxi- 
mate chiral symmetry is spontaneously broken by large 
i/jV- con densates. This phase is denoted by xSB. In this 
regime, quarks have relatively large constituent masses 
which are close to the vacuum values, see Figs. 01 and 
[5J Here the density of all quark flavors is very low and 
even vanishes at T = 0. There is no diquark pairing in 
this phase. The %SB phase is rather insensitive to the 
presence of a nonzero lepton-number chemical potential. 
With increasing fj,L e the phase boundary is only slightly 
shifted to lower values of /z. This is just another mani- 
festation of the strengthening of the 2SC phase due to 
neutrino trapping. 

With increasing temperature, the ^^"Condensates melt 
and the xSB phase turns into the normal quark matter 
(NQ) phase where the quark masses are relatively small. 
Because of the explicit breaking of the chiral symmetry 
by the current quark masses, there is no need for a phase 
transition between the two regimes. In fact, at low chemi- 
cal potentials we find only a smooth cross-over, whereas 
there is a first-order phase transition in a limited region, 
320 < < 360 MeV. In contrast to the xSB regime, the 
high-temperature NQ phase extends to arbitrary large 
values of /i. All these qualitative features are little af- 
fected by the lepton-number chemical potential. 

The third regime is located at relatively low tempera- 
tures but at quark chemical potentials higher than in the 
%SB phase. In this region, the masses of the up and down 
quarks have already dropped to values well below their 
respective chemical potentials while the strange quark 
mass is still large, see left columns of panels in Figs. 01 
and [SI As a consequence, up and down quarks are quite 
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abundant but strange quarks are essentially absent. 

It turns out that the detailed phase structure in this 
region is very sensitive to the lepton-number chemical 
potential. At /j,l b = 0, as discussed in Ref. Q, the pair- 
ing between up and down quarks is strongly hampered 
by the constraints of neutrality and f3 equilibrium. As 
a consequence, there is no pairing at low temperatures, 
T < 10 MeV, and g2SC-type pairing appears at moder- 
ate temperatures, 10 MeV <T<T C with the value of T c 
in the range of several dozen MeV, see the upper panel in 
Fig. [3J The situation changes dramatically with increas- 
ing the value of the lepton-number chemical potential. 
Eventually, the low-temperature region of the normal 
phase of quark matter is replaced by the (g)2SC phase 
(e.g., at {J, = 400 MeV, this happens at fi La ~ 110 MeV). 
With fii, e increasing further, no qualitative changes hap- 
pen in this part of the phase diagram, except that the 
area of the (g)2SC phase expands slightly. 

Finally, the region in the phase diagram at low temper- 
atures and large quark chemical potentials corresponds 
to phases in which the cross-flavor strange-nonstrange 
Cooper pairing becomes possible. In general, as the 
strength of pairing increases with the quark chemical po- 
tential, the system passes through regions of the gap- 
less uSC (guSC), uSC, and gCFL phases and finally 
reaches the CFL phase. (Of course, the intermediate 
phases may not always be realized.) The effect of neu- 
trino trapping, which grows with increasing the lepton- 
number chemical potential, is to push out the location 
of the strange-nonstrange pairing region to larger values 
of ji. Of course, this is in agreement with the general 
arguments in Sec. 11111 

It is interesting to note that the growth of the 
strangeness content with increasing quark chemical po- 
tential could indirectly be deduced from the behavior 
of the electric charge chemical potential hq at T = 0, 
see the solid lines in the right panels in Figs. ^and|3 
The value of /iq reaches its minimum somewhere in a 
range of values of the quark chemical potential around 
/i ~ 440 MeV. This corresponds to the point where the 
strange quark chemical potential /j, s ~ fi — /iQ/3 reaches 
the value of the strange quark mass (see left panels). 
Hence, there are essentially no strange quarks at lower 
values of fi, and a rapidly increasing amount of strange 
quarks at higher values of /i. Since the latter contribute to 
the electric neutralization, this is a natural explanation 
for the drop of \/j.q\ above this point. 

As we mentioned earlier, the presence of the lepton- 
number chemical potential \ij Je leads to a change of 
the quark Fermi momenta. This change in turn affects 
Cooper pairing of quarks, facilitating the appearance of 
some phases and suppressing others. As it turns out, 
there is also another qualitative effect due to a nonzero 
value of HL, e - In particular, we find several new variants 
of gapless phases which did not exist at vanishing 
In Figs. Eland [7| these are denoted by the same names, 
g2SC or gCFL, but with one or two primes added. 

We define the g2SC as the gapless two-flavor color- 



TABLE I: Definitions of gapless phases. 



Name 



Gapless mode(s) 
e(fc) ~ \k-kf\ 



Diquark condensate(s) 



S 2SC 

j2SC' 

juSC 

jCFL 

jCFL' 

?CFL" 



ru-gd, gu-rd 
rd-gu, gd-ru 
rs-bu 
gs-bd 
rs-bu 

gs-bd, rs-bu 



A 3 
A 3 

A 2 , A 3 
Ai, A 2 , A 3 
Ai, A 2 , A 3 
Ai, A 3 , A 3 



superconducting phase in which the gapless excitations 
correspond to rd-gu and gd-ru modes instead of the 
usual ru-gd and gu-rd ones, i.e., u and d flavors are 
exchanged as compared to the usual g2SC phase. The 
g2SC phase becomes possible only when the value of 
(m™ — = (m<2 + A*3)/2 is positive and larger than 

A3. The other phases are defined in a similar manner. 
In particular, the gCFL' phase, already introduced in 
Sec. IIIII is indicated by the gapless rs-bu mode, while 
the gCFL" phase has both gs-bd (as in the gCFL phase) 
and rs-bu gapless modes. Definitions of all gapless phases 
are summarized in Table [I] 



C. Lepton fraction Yl c 

Our numerical results for the lepton fraction Yj, e are 
shown in Fig. El The thick and thin lines correspond to 
two different fixed values of the lepton-number chemi- 
cal potential, fj, Le = 200 MeV and ^ Lc = 400 MeV, 
respectively. For a fixed value of fiL e , we find that the 
lepton fraction changes only slightly with temperature. 
This is concluded from the comparison of the results at 
T = (solid lines), T = 20 MeV (dashed lines), and 
T = 40 MeV (dotted lines) in Fig. © 

As is easy to check, at T = 40 MeV, i.e., when Cooper 
pairing is not so strong, the fj, dependence of Yl c does not 
differ very much from the prediction in the simple two- 
flavor model in Sec. IIIII By saying this, of course, one 
should not undermine the fact that the lepton fraction 
in Fig. 113 has a visible structure in the dependence on fi 
at T = and T = 20 MeV. This indicates that quark 
Cooper pairing plays a nontrivial role in determining the 
value of Yh e - 

Our numerical study shows that it is hard to achieve 
values of the lepton fraction more than about 0.05 in the 
CFL phase. Gapless versions of the CFL phases, on the 
other hand, could accommodate a lepton fraction up to 
about 0.2 or so, provided the quark and lepton-number 
chemical potentials are sufficiently large. 

From Fig. El we can also see that the value of the lepton 
fraction Yj, e « 0.4, i.e., the value expected at the center 
of the protoneutron star right after its creation, requires 
the lepton-number chemical potential /i£ e in the range 
somewhere between 200 MeV and 400 MeV, or slightly 
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higher. The larger the quark chemical potential /i the 
larger /i£, e is needed. Then, in a realistic construction of 
a star, this is likely to result in a noticeable gradient of the 
lepton-number chemical potential at the initial time. This 
gradient may play an important role in the subsequent 
deleptonization due to neutrino diffusion through dense 
matter. 



D. T— jUi e phase diagram 

Now let us explore the phase diagram in the plane 
of temperature and lepton-number chemical potential, 
keeping the quark chemical potential fixed. Two such 
slices of the phase diagram are presented in Fig. The 
upper panel corresponds to a moderate value of the quark 
chemical potential, fi = 400 MeV. This could be loosely 
termed as the "outer stellar core" phase diagram. The 
lower panel in Fig. [7J corresponds to fi = 500 MeV, and 
we could associate it with the "inner stellar core" case. 
Note, however, that the terms "inner core" and "outer 
core" should not be interpreted literally here. The cen- 
tral densities of (proto-) neutron stars are subject to large 
theoretical uncertainties and, thus, are not known very 
well. In the model at hand, the case \i = 400 MeV ("outer 
core") corresponds to a range of densities around 4poj 
while the case /i = 500 MeV ("inner core") corresponds 
to a range of densities around 10po- These values are of 
the same order of magnitude that one typically obtains 
in models (see, e.g., Ref. |3(ilp. 

In addition to the phase diagrams, we also show the 
results for the dynamical quark masses, the gap parame- 
ters, and the charge chemical potentials. These are plot- 
ted as functions of temperature in Fig.|5](/i = 400 MeV, 
"outer core" case) and in Fig. (p = 500 MeV, "in- 
ner core" case), for three different values of the lepton- 
number chemical potential, = MeV (upper panels), 
/il c = 200 MeV (middle panels) and fi Lc = 400 MeV 
(lower panels). 

At first sight, the two diagrams in Fig. [7J look so diffe- 
rent that no obvious connection between them could be 
made. It is natural to ask, therefore, how such a dra- 
matic change could happen with increasing the value 
of the quark chemical potential from /i = 400 MeV to 
fi = 500 MeV. In order to understand this, it is useful 
to place the corresponding slices of the phase diagram in 
the three-dimensional diagram in Fig. [21 

The fi = 500 MeV diagram corresponds to the right- 
hand side surface of the bounding box in Fig. |2 This 
contains almost all complicated phases with strange- 
nonstrange cross-flavor pairing. The fj, = 400 MeV di- 
agram, on the other hand, is obtained by cutting the 
three-dimensional diagram with a plane parallel to the 
bounding surface, but going through the middle of the 
diagram. This part of the diagram is dominated by the 
2SC and NQ phases. Keeping in mind the general struc- 
ture of the three-dimensional phase diagram, it is also not 
difficult to understand how the two diagrams in Fig. [7J 



transform into each other. 

Several comments are in order regarding the zero- 
temperature phase transition from the CFL to gCFL' 
phase, shown by a small black triangle in the phase dia- 
gram at /i = 500 MeV, see the lower panel in Fig. [7J The 
appearance of this transition is in agreement with the 
analytical result in Sec. II I II Moreover, the critical value 
of the lepton-number chemical potential also turns out to 
be very close to the estimate in Eq. (l2*3)l . Indeed, by tak- 
ing into account that M s w 214 MeV and A 2 « 76 MeV, 
we obtain ^ = A 2 - M s 2 /(2p) « 30 MeV which agrees 
well with the numerical value. 

In order to estimate how the critical value of \ij Jc 
changes with decreasing the quark chemical potential be- 
low fi = 500 MeV, we can use the zero-temperature nu- 
merical results for M s and A 2 from Ref. [5|. Then, we 
arrive at the following power-law fit for the ^-dependence 
of the critical value: 

^ 457.4) (622.1 -//), (26) 

for 457.4 < fj, < 500 MeV (in the fit, ^ and n are 
measured in MeV). Note that the CFL phase does not 
appear at T = when fi < 457 .4 MeV @. With in- 
creasing the values of the quark chemical potential above 
/i = 500 MeV, one expects that the critical value of 
should continue to increase for a while, and then de- 
crease when the effects of the cutoff start to suppress 
the size of the gap A 2 . However, the validity of the fit in 
Eq. (|26|l is questionable there because no numerical data 
for /i > 500 MeV was used in its derivation. 

Before concluding this subsection, we should mention 
that a schematic version of the phase diagram in T— /x^ 
plane was earlier presented in Ref. see the left panel 
in Fig. 4 there. In Ref. @, the value of the quark che- 
mical potential was /i — 460 MeV, and therefore a direct 
comparison with our results is not simple. One can see, 
however, that the diagram of Ref. Q fits naturally into 
the three-dimensional diagram in Fig. |3 Also, the dia- 
gram of Ref. @ is topologically close to our version of 
the phase diagram at /j, = 500 MeV, shown in the lower 
panel of Fig. [7] The quantitative difference is not surpris- 
ing: the region of the (g)CFL phase is considerably larger 
at ^ = 500 MeV than at fx = 460 MeV. 



V. CONCLUSIONS 

In this paper, we studied the effect of neutrino trap- 
ping on the phase diagram of neutral three-flavor quark 
matter within the N JL model of Ref. [25| . The results are 
obtained in the mean-field approximation, treating con- 
stituent quark masses as dynamically generated quanti- 
ties. The overall structure of the phase diagram in the 
space of three parameters, namely temperature T, quark 
chemical potential /z, and lepton-number chemical poten- 
tial /iL,., is summarized in Fig. EI This is further detailed 
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in several two-dimensional slices of the phase diagram, in- 
cluding diagrams in the plane of temperature and quark 
chemical potential (see Fig.^J and in the plane of tempe- 
rature and lepton-number chemical potential (see Fig.[7J). 

By making use of simple model-independent argu- 
ments, as well as detailed numerical calculations in the 
framework of an NJL-type model, we find that neutrino 
trapping helps Cooper pairing in the 2SC phase and sup- 
presses the CFL phase. In essence, this is the consequence 
of satisfying the electric neutrality constraint in the quark 
system. In two-flavor quark matter, the (positive) lepton- 
number chemical potential fii, e helps to provide extra 
electrons without inducing a large mismatch between the 
Fermi momenta of up and down quarks. With reducing 
the mismatch, of course, Cooper pairing gets stronger. 
This is in sharp contrast to the situation in the CFL 
phase of quark matter, which is neutral in the absence of 
electrons. Additional electrons due to large can only 
put extra stress on the system. 

In application to protoneutron stars, our findings in 
this paper strongly suggest that the CFL phase is very 
unlikely to appear during the early stage of the stellar 
evolution before the deleptonization is completed. If color 
superconductivity occurs there, the 2SC phase is the best 
candidate for the ground state. In view of this, it might be 
quite natural to suggest that matter inside protoneutron 
stars contains little or no strangeness (just as the cores 



of the progenitor stars) during the early times of their 
evolution. In this connection, it is appropriate to recall 
that neutrino trapping also suppresses the appearance of 
strangeness in the form of hyperonic matter and kaon 
condensation [3(J • Our results, therefore, are a special 
case of a generic property. 

After the deleptonization occurs, it is possible that the 
ground state of matter at high density in the central re- 
gion of the star is the CFL phase. This phase contains a 
large number of strange quarks. Therefore, an abundant 
production of strangeness should happen right after the 
deleptonization in the protoneutron star. If realized in 
nature, in principle this scenario may have observational 
signatures. 
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FIG. 2: General structure of the phase diagram of neutral dense quark matter in the three-dimensional space spanned by the 
quark chemical potential fi, the lepton-number chemical potential [iL e , and the temperature T. 
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FIG. 3: The phase diagrams of neutral quark matter at fixed lepton-number chemical potentials = (upper panel), 
(iLe — 200 MeV (middle panel), and = 400 MeV (lower panel). Note that the upper diagram, describing matter in absence 
of neutrino trapping, is the same as in Fig. 1 of Ref. Q. 
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FIG. 4: Dependence of the quark masses, of the gap parameters, and of the electric and color charge chemical potentials on the 
quark chemical potential at a fixed temperature, T = MeV (three upper panels) and T = 40 MeV (three lower panels). The 
lepton- number chemical potential is kept fixed at fj,L e = 200 MeV. 
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FIG. 5: Dependence of the quark masses, of the gap parameters, and of the electric and color charge chemical potentials on the 
quark chemical potential at a fixed temperature, T = MeV (three upper panels), and T = 40 MeV (three lower panels). The 
lepton- number chemical potential is kept fixed at [iL a = 400 MeV. 




FIG. 6: Dependence of the electron family lepton fraction Yl c for /j,L e — 200 MeV (thick lines) and /iL e = 400 MeV (thin lines) 
on the quark chemical potential at a fixed temperature, T — MeV (solid lines), T = 20 MeV (dashed lines), and T = 40 MeV 
(dotted lines). 
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FIG. 7: The phase diagrams of neutral quark matter in the plane of temperature and lepton-number chemical potential at two 
fixed values of quark chemical potential: /i = 400 MeV (upper panel) and /i = 500 MeV (lower panel). The triangle in the lower 
panel denotes the transition point from the CFL phase to the gCFL' phase at T = 0. 
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FIG. 9: Dependence of the quark masses, of the gap parameters, and of the electric and color charge chemical potentials on 
the temperature at a fixed value of the lepton-number chemical potential, Hl c ~ MeV (three upper panels), fJ,L e = 200 MeV 
(three middle panels), and fj,L e = 400 MeV (three lower panels). The quark chemical potential is fj, = 500 MeV which may 
correspond to the conditions at the center of a quark core inside a neutron star. 



